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ABSTRACT 

SCUBA-2 is a 10000-bolometer submillimetre camera on the James Clerk Maxwell Tele- 
scope (JCMT). The instrument commissioning was completed in September 2011, and full 
science operations began in October 2011. To harness the full potential of this powerful new 
astronomical tool, the instrument calibration must be accurate and well understood. To this 
end, the algorithms for calculating the line-of-sight opacity have been improved, and the de- 
rived atmospheric extinction relationships at both wavebands of the SCUBA-2 instrument are 
presented. The results from over 500 primary and secondary calibrator observations have al- 
lowed accurate determination of the flux conversion factors (FCF) for the 850 and 450 jam 
arrays. Descriptions of the instrument beam-shape and photometry methods are presented. 
The calibration factors are well determined, with relative calibration accuracy better than 5 
per cent at 850 jam and 10 per cent at 450 jam, reflecting the success of the derived opacity 
relations as well as the stability of the performance of the instrument over several months. The 
sample-size of the calibration observations and accurate FCFs have allowed the determination 
of the 850 and 450 \im fluxes of several well-known submillimetre sources, and these results 
are compared with previous measurements from SCUBA. 

Key words: submillimetre - atmospheric effects. 



1 INTRODUCTION 

Ground-based submillimetre astronomy has progressed markedly 
in the last fifteen years. SCUBA-2 ( [Holland et aLl|20T3]) is the 
successor to the highly successful SCUBA instrument (* Holland| 
[eTal. 1999 ) which was operational on the JCMT from 1998 to 
2005. In both cases, a limitation of ground-based submillimetre 
observing lies in the attenuation and distortion of the incoming 
radiation by the atmosphere. The dominant source of this attenua- 
tion is absorption by water vapour which, even on high sites such 
as Mauna Kea, limits submillimetre observations to a few narrow 
wavebands between the strongest absorption lines. 

The opacity in these bands can vary rapidly and markedly 
as a function of precipitable water vapour (PWV), particularly 
at shorter wavelengths. Rapid measurements of the line-of-sight 
PWV are required in order to compensate adequately for this 
attenuation in the observation. A model of the atmospheric 
transmission as a function of wavelength must be used to relate 
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the PWV to the opacity in the filter bands of the camera. This 
relation must be determined accurately. Uncertainties affect the 
short-wavelength windows significantly as they are more opaque 
and are more highly attenuated as the weather becomes poorer. 
These atmospheric effects have been investigated previously for 
SCUBA ( [Archibald eraLl|2002t and for Bolocam ( [Sayers et alT] 
|20T0l ). 

The second parameter required to calibrate submillimetre 
continuum data is a measure of the optical throughput of both the 
telescope and instrument itself. This value is referred to as the 
flux conversion factor (FCF) which converts the measured power 
in picowatts (pW) into astronomical fluxes in janskys (Jy). To 
obtain this factor, a series of astronomical sources of known flux 
are observed, preferably in a range of atmospheric conditions and 
over an extended time. Variations in this factor are dominated by 
optical effects from thermal distortions of the dish (particularly in 
the early evenings). Poor focus or even submillimetre seeing in 
the atmosphere will also result in deviations in the peak flux of a 
compact source. 

This paper describes the atmospheric transmission relations 
for both SCUBA-2 filter bands and how these are applied to 
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Figure 1. Examples of the ~3 pW peak-to-peak fast ramps of the heater 
power that are used to calculate the flat-field response (right) of each in- 
dividual TES bolometer. The heater is ramped eight times in a ten second 
interval. The top plots show a typical, well behaved bolometer ramp and 
the corresponding current change as a function of heater power, with the 
responsivity determined from a fit to the inverted slope dl/dP. Bolometers 
are flagged and removed as bad if showing too poor a S/N (middle plots) or 
a severe non-linear response (bottom). 

the data. The calibration sources are listed and the methods of 
observation and data reduction are described. The results of aper- 
ture photometry used to calculate the flux conversion factors are 
presented. The beam shapes at both wavelengths are modeled and 
discussed. Finally, a number of well-known submillimetre sources 
are presented with their calibrated fluxes at 850 ^tm and 450 |^m. 
These values are compared with previous results measured with 
SCUBA. 



2 SCUBA-2 

The commissioning of SCUBA-2 was completed on the JCMT in 
September 201 1. The camera observes in two wavebands, 850 and 
450 |j,m simultaneously, with 5120 pixels in each focal plane. The 
full description of the instrument and its performance characteris- 
tics are described in [Holland et al.| ( [20T 3 ). Results from the instru- 
ment characterisation and on-sky commissioning are presented in 
[Bintley et al.| ( [20T0|[20T2l ) and |Dempsey etaLl ( [20T2l ). The SCUBA- 
2 bolometers are integrated arrays of superconducting transition- 
edge sensors (TESs). The devices are voltage biased to a tempera- 
ture near the centre of their transition range, whereupon any tem- 
perature fluctuation results in a current change through the TES. 
The resulting magnetic field change is then amplified using chains 
of superconducting quantum interference devices (SQUIDs) and 
the amplified current is digitised. SCUBA-2 detectors have indi- 
vidually coupled resistive heaters which are used to compensate 
for changes in incident optical power, thereby allowing the detec- 
tor bias point to remain constant over a wide-range of sky powers. 

2.1 Flat-fielding 

On-sky observations with SCUBA-2 proceed following an array 
setup, which is conducted with the internal instrument shutter 
closed to optimise the SQUID settings. When the shutter is opened 
the power from the resistive heater coupled to each bolometer 
is decreased to compensate for the incident optical power. This 
adjustment is termed a "heater track". The resistor coupling to 
the TES varies between sub-arrays. The "effective" resistance 
can be calculated by measurement of the I-V curve obtained at 



different heater power settings. These resistances were shown to 
be significantly different for each individual sub-array. The relative 
optical response of the arrays was determined by measurement 
of the signal when observing bright sources such as Jupiter. The 
coupling factors were adjusted such that the picowatt response 
to the planet signal was equal in all sub-arrays. The I-V curves 
were used to obtain the average absolute pW calibration of the 
bolometers and then the flux-determined factors were used to 
calibrate the relative sub-array response, as this is ultimately the 
measure of the optical response through the entire system from an 
astronomical source. 

The heaters are then employed to calibrate the individual 
bolometer responsivity. The response of the detectors is deter- 
mined by measuring the slope of the current/power change as the 
heater power below each pixel is ramped quickly over a small 
range (a few pW peak-to-peak), as shown in the left-hand plots 
in Figure [T] The flat-field response is the inverse gradient of the 
current change as a function of heater power (or dl/dP) of each 
individual bolometer. At this point the quality of the bolometer 
performance is assessed: bolometers with responsivities above or 
below a threshold limit are removed (these limits are extreme, and 
are applied to remove the handful of bolometers which display 
responsivities too high or low to be physically accurate), along 
with bolometers whose response, dl/dP, is non-linear, or shows 
too poor a S/N. Figure [T] shows a good bolometer trace (top) and 
example of a bolometer that was flagged resulting from a poor S/N 
(centre) and extreme non-linear response (bottom). A flat-field is 
performed at the start of every science observation, and the quality 
flags typically remove less than 10 per cent of the total bolometer 
yield at 450 ^tm and less than 20 per cent of the array yield at 
850 i^m with the percentage and pattern of the subsequent masking 
being typically consistent. 

Initially, this flat-fielding was completed in the dark, with the 
internal instrument shutter closed, assuming that the heater track 
as the shutter opened would compensate sufficiently for the change 
in optical power. Testing showed, however, that better bolometer 
yields and more stable operation was achieved if this flat-fielding 
was conducted again once the shutter was open, particularly when 
there was sky instability. Dark flat-fields are still conducted at the 
end of an array setup to check if the setup was successful. 



3 SUBMILLIMETRE ATMOSPHERIC TRANSMISSION 

The atmosphere limits ground based observations at submillimetre 
wavelengths, providing only a few semi-transparent windows even 
at a high, dry site such as the summit of Mauna Kea. SCUBA-2, like 
its predecessor SCUBA, has been designed to take advantage of the 
850 and 450 \im atmospheric windows. The measured SCUBA-2 
filter profile at both wavelengths are shown superimposed on the 
submillimetre atmospheric transmission windows in Figure |2] The 
filter profiles were measured using a Fourier transform spectrome- 
ter on a single bolometer at a time. The filter profiles were measured 
on a number of bolometers moving away from the centre of the fo- 
cal plane and no discernible difference was observed in the profile. 

Atmospheric attenuation at a given waveband in the submil- 
limetre can be described as follows, assuming a plane-parallel at- 
mosphere: 
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Figure 2. 450 fxm (blue solid) and 850 \im (red dotted) SCUBA-2 filter 
profiles superimposed on the submillimetre atmospheric transmission curve 
for Mauna Kea assuming 0.5 mm of precipitable water vapour. 
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where /m and /o are the measured signal at the telescope and at the 
top of the atmosphere respectively, t is the extinction coefficient, 
and A is the airmass (eg. [Stevens & Robson| ( |199 4 )). To accurately 
reconstruct the un-attenuated signal /o, the extinction coefficient x 
must be well described. At the shorter of SCUBA-2's filter bands 
a 20 percent error in t can result in a corresponding flux error of 
50-80 percent ( [Archibald et al.|2002| ). 
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Figure 3. CSO 225 GHz radiometer data vs precipitable water vapour (in 
millimetres) as derived by the 183 GHz line-of-sight water vapour monitor 
at JCMT. WVM and CSO data collected in the most stable part of the night, 
between 9 pm and 3 am, is shown. The full range of PWV fitted is shown on 
the left, the red rectangle highlighting the best weather conditions (< 3 mm 
PWV) which is shown in detail in the right hand plot. 

( [Dempsey & Friberg|200"8l ). 

If we adjust the PWV to its zenith value, PWVzen by dividing 
by the airmass, there exists a simple linear relation between it and 
the opacity at 225 GHz. We can then directly compare the opacity 
derived in this way by the WVM to the independently produced 
CSO 225 GHz opacity. The WVM and CSO values show excellent 
correlation for the stable part of the night between 9 pm and 3 am, 
with deviations occurring in the early evening and at sunrise when 
conditions become more variable. The PWVzen is related to the 
opacity at 225 GHz, t(225)_^ by : 



T(225), 



0.04 PW^Vzen + 0.017. 



(2) 



where PWVzen is in units of millimetres. The comparison 
between the CSO 225 GHz opacity at zenith and the PWVzen from 
the WVM for a full year of observing in 201 1 is shown in Figure[3] 



3.1 The Water Vapour Monitor 

The extinction coefficient, T, must be explicitly determined for 
each observing waveband. The Caltech Submillimetre Observatory 
(CSO) has a fixed-azimuth 225 GHz tipping radiometer that has 
been in operation throughout the lifetime of both SCUBA and 
SCUBA-2 ( [Radford|20Tlj ). This completes a skydip every fifteen 
minutes and calculates a zenith opacity using a plane-parallel 
atmospheric model. The wavelength dependent extinction re- 
lations for SCUBA were calculated as a function of the CSO 
tipper 225 GHz opacity ( [Archibald et al.|2002| ). The JCMT has a 
183 GHz water vapour monitor (WVM) installed in the receiver 
cabin that measures the precipitable water vapour (PWV) along 
the line-of-sight at 1.2-second intervals { [Wiedner et al.[|200T| ). 
The WVM has the advantage of better time resolution, which is 
important given that SCUBA-2 collects data at a rate of 200 Hz. 
Also, by measuring directly along the line-of-sight of the telescope 
we are not limited to the fixed azimuth of the CSO tipper, and 
we do not need to assume a plane-parallel atmosphere. Skydips 
taken with the 345 GHz HARP instrument ( [Buckle et ar]|200^ 
on the JCMT, as well as previous comparisons between SCUBA 
skydips and the CSO tipper show that the atmosphere deviates 
from a plane-parallel approximation about 30 per cent of the time, 
particularly in the early evening when the atmosphere is unstable 



3.2 Tau relations 

The extinction coefficients at each SCUBA-2 wavelength will be 
described in terms of the PWV determined by the water vapour 
monitor. The relation was derived by collating the calibration ob- 
servations of sources of known flux taken over the six month period 
from 2011 May until 2011 October, and then in additional science 
time until 2012 May. These observations were reduced without ap- 
plying any assumed extinction correction and then fitted using a 
least- squares algorithm to the following formula, where the precip- 
itable water vapour (PWV) is measured in millimetres: 

TA = a(PWV-6). (3) 

The resulting extinction coefficients at each wavelength are: 

T850 = 0.179 (PWV + 0.337) (4) 

and 

T450 = 1.014 (PWV + 0.142). (5) 

The data from which these extinction coefficients were deter- 
mined is presented and discussed in Section [54]below. The quality 
of these fits is thus determined by the scatter and systematics in the 
calibrator data, and is accounted for in the resulting error quoted 
for the flux conversion factor in Section [54l The SCUBA-2 extinc- 
tion correction coefficients in terms of the opacity at 225 GHz are: 
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Figure 4. The top left image shows the inner 120 arcseconds of the co- 
added map of Uranus at 850 |j.m, with contours in square-root scaling from 
the peak down to 5 per cent of the peak flux. The top right image shows the 
inner 60 arcseconds of the corresponding Uranus map at 450 |Lim. The lower 
plots show the low-level, large scale beam pattern at both wavelengths, with 
the colour scaled from 0. 1 per cent of the peak flux down to zero at both 
wavelengths. 



T850 = 4.6 (t225 - 0.0043) andT45o = 26.0 (T225 - 0.012), to allow 
easy comparison with the previous extinction terms for SCUBA. 



4 SCUBA-2 BEAM SHAPE 

The JCMT beam as seen with SCUBA-2 has been determined from 
good quality maps of Uranus obtained between 2011 May and 
2012 June. The maps were reduced using the method described 
later in Section |5.2| The maps co-add the measured signal from 
all working bolometers. Individual bolometer scans of Uranus 
show that the beam profile varies from the centre to the edge of 
the focal-plane by less than 5 per cent. The individual maps were 
then co-added after adjustment for pointing differences. The size 
of Uranus was ^3.5arcsec during this period. This reduction 
filters map structures on scales approximately larger than 5 arcmin 
( [Chapin et al. ( 2013 )). Thus any beam contribution larger than this 
will not be seen in SCUBA-2 maps. 

The beam is described here in terms of two Gaussian compo- 
nents - a main-beam Gaussian, Gmb, and an additional error-beam 
component, Gs describing the response due to large scale surface 
distortion and diffraction (Airy pattern). While the main beam is 
well described by a Gaussian this does not apply to the second 
component. However, using a Gaussian is a simple way to capture 
the volume of this error-beam. A more complete description of 
the second component would by asymmetric and stationary in 
azimuth and elevation as a result of the mounting of the telescope. 
It is also well known that the point source efficiency is affected by 
thermal surface distortions early in the evening and mid-morning. 
The associated aberrations would increase and possibly extend the 
second component. Drops of aperture efficiency as large as 50 per 
cent have been reported in the 450 |^m heterodyne systems, and 



Figure 5. The upper plots show the azimuth (cyan) and elevation (blue) 
cuts through the centre of the Uranus beams from Figure |4] with the two- 
component fitted model overlaid in red. The lower panels show a zoom in 
to the secondary component of the beam. The slice through the 850 |Ltm in 
particular, shows distinct asymmetry that is not well approximated by the 
Gaussian model. 



SCUBA measured FCF variations of up to 30 per cent at 850 \im 
under such conditions ( [Jenness et al.|2002) . 

Figure |4] shows the co-added Uranus maps at both wave- 
lengths. The upper plots show the inner detail of the beam pattern, 
while the lower plots show the low-level emission at larger scales. 
Normalised radial slices through the beam in the azimuth (blue) 
and elevation (cyan) directions are shown in Figure |5] More dis- 
tinct asymmetry is noted in the 850 \xm cuts than at 450 |^m. The 
beam shape was then fitted as the normalised sum of the two Gaus- 
sian components as given by: 



Gt 



= a Gmb + /3 Gs 



(6) 



where the Gaussian profile, G, is described as exp[— 4 In 2(r/9)^] 
where r is the radial distance from the centre in arcseconds and 9 
is the full- width half maximum (FWHM) in arcsec. The result of 
the Gaussian fitting is given in Table [T] The integral of the profile 
provides the volume. 



(7) 



4/n(2) 

in arcsec^. The widening effect resulting from the finite size of 
Uranus has been removed from the widths given in the table. The 
main beam widths are close to theoretical values estimated with 
a 3 dB taper - these are 7.4 and 12.75 arcsec, respectively. This 
includes widening caused by the finite pixel size F/2A at 850 ^im 
and F/Aat450M,m. 

Apart from not estimating the taper correctly the measured 
beam is also affected by the atmosphere (seeing), telescope 
focusing errors as well as large scale surface errors. In addition 
slight pointing shifts between the co-added maps will also widen 
the beam even if the maps are aligned with each other before 
co-adding. Despite these factors, the volume of the error-beam is 
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Table 1. Results of two component Gaussian fits 
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satisfactorily described by this model. Subtracting the modeled 
main-beam Gaussian from the co- added Uranus map in Figure [4] 
and integrating over the remaining volume allows empirical com- 
parison of the error-beam with that described by the model. The 
measured error-beam is 24 and 39 per cent of the total integrated 
flux at 850 and 450 \im, respectively, compared to 25 and 40 per 
cent predicted by the model. The size of the second Gaussian 
is rather close to the first maximum in the Airy pattern from a 
uniform aperture. The maxima are roughly at diameters 40 and 20 
arcseconds for 850 and 450 ^tm, respectively. 

Strong point sources observations clearly show a low-level 
emission ring at large scales (see Figure |4]). The feature results 
from a slight mis-match between the telescope focal length and 
the focal length of the individual panels. The ring is 350 and 183^^ 
in diameter at 850 and 450 ^tm, respectively. Though visible, the 
contribution to the beam integral is less than 1 per cent (850 \im) 
and 3 percent (450 ^im). This contribution is not included in 
Table[T] 



5 FLUX CALIBRATION 

Nightly flux calibration is achieved by observation of astronomical 
sources with known flux properties. To convert the measured signal 
of a source from picowatts to janskys a flux conversion factor must 
be applied. The primary factors affecting the FCF are the optical 
path to the detectors, including the shape of the dish, and the filter 
profiles at each wavelength. If the optical quality is stable, and flat- 
fielding of the instrument and extinction correction are properly 
calculated, the FCF should be a constant factor. 
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Figure 6. Aperture photometry curve of growth normalised for a 60 arcsec 
aperture, at both wavelengths. 



additional potential secondary calibrators were also observed and 
are listed with their SCUBA fluxes in Table [2] 

5.2 Observation method and data reduction 

All calibrators were observed using a daisy scan-pattern ( |Holland| 
|et al.|[20T3] ), designed specifically for observation of point and 
compact sources, with observations lasting approximately four 
minutes each. The data were reduced using the SMURF routine 
MAKEMAP ( [Chapin et"ar||2013| [Jenness et~aL]|2011| ), using a 
parameter file designed specifically for bright point sources. 
The reduction of compact calibrators uses a combination of 
common-mode rejection and high-pass filtering of the bolometer 
data (which effectively suppresses information on scales larger 
than 3.3 arcmin) to remove low-frequency noise. The ringing one 
might expect from such filtering is then controlled by constraining 
the map to a value of zero beyond 60 arcsec from the source (well 
beyond the bulk of the power in the central lobe) for all but the 
final iteration. For point-like sources, the gain (independent of 
the brightness of the source) is effectively one. These parameters 
are discussed in detail in Section 4.1 of [Chapin et al.|pQ131 l. The 
calibrations were taken between 2011 May and 2012 May, and 
over 500 observations (at each wavelength) were included in this 
analysis. Maps were regridded to a 1 -arcsec pixel scale, to optimise 
the quality of the fits to the peak of the beam. 



5.1 Calibration sources 

Mars and Uranus were the primary calibrators. The brightness tem- 
peratures of these planets were obtained using the JCMT FLUXES 
software ( Privett et al. 2005) based on models from Wright|((T976| 
and Moreno (2010|^respectively. It should be noted that the newer 
Uranus brightness temperatures differ from those used for SCUBA 
(a decrease of approximately 5 per cent and 2.5 per cent at 850 ^im 
and 450 |^m respectively). | Wright] (| 1976)) estimates an uncertainty 
in the Mars model flux of ±5 per cent. The Uranus model is quoted 
as having an absolute uncertainty of 5 per cent. Secondary calibra- 
tors used on a nightly basis were CRL 618 and CRL 2688. A set of 

^ Moreno, R. "Neptune and Uranus planetary brightness temperature 
tabulation" ESA Herschel Science Centre, ftp://ftp.sciops.esa.int/pub/hsc- 
|calibration[ 2010 



5.3 Aperture photometry 

The total integrated flux for each calibrator map was calculated in 
a 60 arcsec diameter aperture, with the background level calculated 
within an annulus of inner diameter 90 arcsec and outer diameter 
120 arcsec (1.5 and 2.0 times the source aperture). As discussed 
in Section |4] the SCUBA-2 beam has a broad error-beam, and a 
60 arcsec aperture does not fully encompass the beam, particularly 
at 850 i^m. However, increasing the aperture, and thus by necessity, 
the annulus, beyond this size showed a significant increase in 
noise, and produced larger scatter in the results. 

The 60 arcsec aperture was chosen as a compromise between 
minimising the scatter in the results whilst optimising the total 
flux measurement. The curve of growth in Figure |6] shows the 
ratio of Aperture flux to Total Flux for a co-add of several Uranus 
observations, normalised to an aperture of 60 arcsec. At 450 \im. 
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Table 2. Previously published flux densities of SCUBA-2 potential secondary calibrators. 



SOURCE 


RA 


Dec. 


850 [I] 


450 [I] 


850 [P] 


450 [P] 


ref 




(J2000) 


(J2000) 


(Jy) 


(Jy) 


(Jy beam"-*^) 


(Jy beam" ^ ) 




CRL618 


04:42:53.67 


+36:06:53.17 


4.73 ± 0.33 


12.1 ±2.2 


4.55 ± 0.2 


11.5 ± 1.5 


d[^],b[^] 


CRL2688 


21:02:18.27 


+36:41:37.00 


6.39 ±0.51 


30.9 ± 3.8 


5.9 ± 0.2 


24.0 ±2.1 


d[^],b[^] 


Arp220 


15:34:57.27 


+23:30:10.48 


0.83 ± 0.086 




0.7* 


5.0* 


f[^],c[^] 


PV Cep 


20:45:54.39 


+67:57:38.8 


1.0 ±0.02 


6.5 ± 0.015 






a 


MWC 349 


20:32:45.53 


+40:39:36.6 


2.6 ± 0.007 


5± 1.1 






a 


HD169142 


18:24:29.78 


-29:46:49.37 


0.565 ± 0.01 


3.34 ±0.115 






a 


V883 Ori 


05:38:19 


-07:02:2 


1.41 ±0.022 


9.45 ±0.17 








HLTau 


04:31:38.44 


+18:13:57.65 


2.36 ± 0.24 


9.9 ± 2.0 


2.3 ±0.15 


10.0 ± 1.2 


d[^]^b[^] 


BVP 1 


17:43:10.37 


-29:51:44.00 






1.069 ±0.02 


11.318 ±0.6 


e 



[I] indicates integrated flux measurements. [P] indicates peak flu x measurem ents. References: (a) Sandell et al. ( 2011 », background-corrected 
integrated flux, (b) SandeU (2003 1, peak flux in Jy beam"!, (c) Truch et al.'(2008 i, BLAST primary calibrator. *Peak flux based on fit to SED 
curve fitted to measurements at other wavelengths. (d) Jenness et al. (2002) integrated flux in a 40 arcsec aperture, (e) [Barnard et al.H2004) 
Peak flux. Denoted as object G6 in that text. (f) |Lisenfeld et al.| ( |2000y integrated flux. (g) |Sandell & Weintraub| ( |200l} , background-corrected 
integrated flux. 



Table 3. When integrating over an aperture of diameter, d, in arcseconds, 
the FCF given in Section [5!4] must be divided by the following factor, given 
in columns denoted 450i? and 850i?. The factors are derived from the mea- 
surements plotted in Figure [6] These factors are only valid for unresolved 
sources. 



d (arcsec) 450^ 850^ 



20 


0.72 


0.69 


25 


0.79 


0.79 


30 


0.85 


0.85 


35 


0.88 


0.88 


40 


0.91 


0.91 


45 


0.94 


0.94 


50 


0.97 


0.97 


55 


0.98 


0.99 


60 


1.00 


1.00 


65 


1.01 


1.01 


70 


1.02 


1.02 


75 


1.02 


1.03 


80 


1.02 


1.04 


85 


1.03 


1.05 


90 


1.03 


1.05 


95 


1.03 


1.06 


100 


1.03 


1.07 


105-120 


1.03 


1.08 


>120 




1.09 



this aperture underestimates the total flux by approximately 4 per 
cent while at 850 \xm this value is approximately 8 per cent. In 
the results presented here, the Uranus calibration maps and the 
secondary sources are both reduced using a 60 arcsec aperture, 
therefore no correction is required in the FCF values quoted below. 
However, if a different aperture size is chosen for a science source, 
the result needs to be divided by the appropriate factor derived 
from Figure |6] and listed for reference in the lookup Table [3] These 
factors are only applicable to photometry of unresolved sources. 

The FCF is simply described as the ratio between the known 
flux S, in Jy, and the measured signal in pW, /o, summed over the 
source. To apply the FCF when calculating a flux integrated over a 
source, the size of the pixels in the output map must also be taken 
into account. This factor is called the FCFarcsec, and is given by: 



S 

FCF arcsec — t 4 i (8) 
io A 

where A is the pixel area in arcsec^, giving the units of 
Jy pW~^ arcsec"^. 

It is sometimes preferable to measure the peak flux of a source, 
rather than the integrated value. To fit the peak of each calibration 
map, a Gaussian-like 2D fit was applied to each calibrator using the 
KAPPA application BEAMFIT ( ,Currie et al.|200 8), with the power 
index allowed to be a variable in the fit. Unlike the total flux cal- 
culated in a large aperture, the source peak is highly susceptible 
to errors in telescope focus, as well as temperature-dependent dis- 
tortions in the dish shape, most noticeable early in evening obser- 
vations. The peak FCF, FCFpeak, is simply the ratio of the known 
peak flux in Jy to the measured peak flux in pW. 

5.4 Results 

Figures [7] [8] and [9] show the results of the FCF calculations for 
the standard calibrators, Uranus, Mars, CRL618 and CRL2688. 
Figure [T] shows FCFarcsec as a function of time (top), and versus 
atmospheric transmission (bottom). At both wavelengths the 
results show no systematic deviation, indicating the instrument and 
optical performance were consistent during this year-long period. 

Figure |8] shows the same comparisons, but for the FCFpeak- 
These results shows larger deviation from the mean, particularly 
at 450 i^m. This is to be expected given this measure is more sus- 
ceptible to focus and temperature variations in the optics, however, 
again, no systematic deviation is observed. Figure [9] shows the his- 
togram distribution of FCFarcsec and FCFpeak for the combined 
set of calibrations at 850 i^m (top left panel) and 450 ^tm (top right 
panel). The resulting mean FCFs at each wavelength are thus: 

FCF [850] ^^^3^^ = 2.34 d= 0.08 Jy pW"' arcsec-' (9) 

and 

FCF [450]^^^^^^ = 4.71 ± 0.5 Jy pW"^ arcsec"^ (10) 

The corresponding peak FCF results are: 
FCF[850]p^^k = 537 ± 26 Jy pW"' (11) 

and 
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Figure 7. Flux conversion factor, FCFarcsec results, at 850 \im (left) and 
450 |j.m (right), plotted as a function of time (top) between May 2011 and 
May 2012, and as a function of atmospheric transmission (bottom). The 
solid line indicates the average value, and the dotted lines denote the 1-cr 
error. 
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Figure 9. The top left plot shows the histogram distribution of the 
FCFarcsec (black, soHd) and FCFpeak (red, outhne) at 850 |Ltm for the data 
presented in Figure^ while the top right plot shows the corresponding 
distributions at 450 |Ltm. The two lower panels show the FCFarcsec versus 
FCFpeak at both wavelengths. A linear fit to this relation allows an empiri- 
cal estimate of the beam area. 
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(12) 



The error is the 1-cr deviation from the mean. The relative flux 
calibration is thus better than 5 per cent at 850 |j,m and ^10 per 
cent at 450 |j,m. In comparison, the relative flux calibration of 
SCUBA from [Jenness et aL| ( [20021 ) was 10 per cent at 850 \im 
and 20 per cent at 450 \im. Current instruments with comparable 
wavebands include LAB OCA, on the APEX telescope (Sir ingoj 
let al. 12009 ) which measures a calibration accuracy of 10 per cent 
at 870 M^m and BOLOCAM, at the Caltech Submillimeter Obser- 
vatory, with an 8 per cent uncertainty at 1.1 millimetres ( jAguirrel 
let al.|2011 |. SHARC-II, the 350 ^im bolometer camera at Caltech 
Submillimeter Observatory, reports a calibration uncertainty of 
15 per cent (Kovacs et al._2006 ), and as it shares the same site as 
SCUBA-2 it provides a point of comparison for the SCUBA-2 
450 ^Lm accuracy, as the 350 ^tm and 450 ^im atmospheric win- 
dows are similarly affected by small variations in PWV. This 
improvement by a factor of two in the calibration accuracy at both 
wavelengths, in comparison to SCUBA, can be attributed to the 
increased time-resolution measurements of the atmospheric condi- 
tions by the WVM, and the stability of the instrument performance. 



Figure 8. Flux conversion factor, FCFpeak results, at 850 |xm (left) and 
450 |xm (right), plotted as a function of time (top) between May 2011 and 
May 2012, and as a function of atmospheric transmission (bottom). The 
solid line indicates the average value, and the dotted lines denote the l-a 
error. 



The absolute calibration accuracy must also take into account 
the uncertainty in the modeled fluxes of our primary calibrators. 
Mars and Uranus, which as noted in Section [5?7] are both ± 5 per 
cent. Less than a percent difference is observed in the mean FCF 
derived for the calibrators using only one or other of these sources 
implying that the estimate of the uncertainty on these models is 
an upper limit and we have a better calibration than the addition 
of this error estimate would suggest. Taking the case that we have 
an overall model error of 5 per cent, the quadrature addition of the 
model error and relative FCF error would give an upper limit on 
the absolute calibration uncertainty of 8 per cent at 850 ^im and 12 
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Figure 10. 850 |xm (left) and 450 |xm (right) image of CRL 618. The colour 
scales linearly from the background level to the peak as indicated by the 
colour bar on the right-hand side of each plot. Contours indicate 1, 2, 5, 10, 
20 and 50 per cent of the peak flux. The 3-a level at 850 |xm is 0.3 per cent 
of the peak flux at 17 mJy, while at 450 fxm 3-a is 1 per cent of the peak at 
116mJy. 



Figure 11. 850 |xm (left) and 450 ^tm (right) image of CRL 2688. The 
colour scales linearly from the background level to the peak as indicated 
by the colour bar on the right-hand side of each plot. Contours indicate 1, 
2, 5, 10, 20 and 50 per cent of the peak flux. The 3-cr level at 850 [im is 0.2 
per cent of the peak flux at 1 1 mJy, while at 450 |Ltm 3-cr is 0.5 per cent of 
the peak at 83 mJy. 



per cent at 450 ^tm. 

The lower panels in Figure |9] compare the FCFpeak versus 
FCFarcsec at cach wavelength. This allows an empirical estimation 
of the beam area, Ae, by fitting the linear relation between the peak 
and integrated values. The beam areas derived in this way therefore 
determine the effective FWHM at each wavelength. Assuming 
a gaussian beam, the effective FWHM, 9e, at each wavelength 
can be determined as y^Ae/1.133), resulting in 9e of 14.1 and 
9.6 arcsec at 850 and 450 \im respectively. Note that this is wider 
than the fit to the main beam component at each wavelength as 
shown in Table [T] as it accounts for the extended emission in the 
shoulders of the beam. The effective FWHM should relate to the 
derived, two-component beam as 9e = \/a(0MB)^ + /3(9s)^. 
The 9e predicted from the two-component model gives 14.6 
and 9.8 arcsec at 850 and 450 |^m respectively, which compares 
reasonably with the estimate derived from the calibrations. 



6 CALIBRATION SOURCE FLUXES 

Several well- studied submillimetre point sources were also regu- 
larly observed as potential new calibrators. Using the new extinc- 
tion corrections and FCFs derived above, these sources were cal- 
ibrated and compared with previous determinations of their flux. 
Table [4] shows the calibrator fluxes for all potential calibrators that 
were observed. In all cases, a co-added map of a large set (or in 
some cases, all) of the observations was produced and the inte- 
grated fluxes were calculated in a 60 arcsec aperture as described 
in Section [53] The table reports the integrated fluxes and also the 
peak fluxes as derived from a single Gaussian fit to the source pro- 
file. 



6.1 CRL 618 

CRL 618 is also extremely well-observed in the submillimetre 
( |Sandell| [20031 [Jenness et al1[20T0l [2002] [Siringo et all[2009l ) 
and is compact, with a measured FWHM of 13.7 and 8.0 arcsec 
at 850 and 450 ^im, respectively. Figure [To| shows CRL 618 at 
both wavelengths. The calculated integrated flux at 850 ^im is 
slightly higher at 5.0 Jy in comparison to results presented by 
[Jenness et aL| p002| ), though the result does fall within their 



error range. The possibility that this source varies over time has 
been suggested ( [Knapp et al.[ |1993). No source variability was 
seen in the SCUBA-2 results between 2011 July and 2012 May, 
which agrees with the results from SCUBA between 1997 and 
2005 as presented by [Jenness eraLlpOTO] ). At 870 ^im LABOCA 
measure a peak flux of 4.9=b0.2Jy beam" ^( Siringo et al.[[2009| ) 
which agrees well with our result of 4.89=b0.24Jy beam" . At 
450 i^m, the integrated and peak fluxes agree extremely well with 
the published SCUBA counterparts. 



6.2 CRL 2688 

CRL 2688 is a well-known compact submillimetre source that has 
been used as a calibrator for SCUBA (Sandell 2003|), and other 
submillimetre instruments including BLAST (Truch et al.[[2008[ l. 
Over 170 observations of the source were taken with SCUBA-2 
during 201 1 and early 2012, and the co-added map of 55 of the best 
weather observations is shown at both wavelengths in Figure[TT] At 
850 i^m, the source is compact, though possibly slightly resolved, 
and the measured FWHM is 14.2 arcsec. At 450 ^tm, we measure a 
FWHM of 9.4 arcsec indicating the source is extended. This is in 
agreement with previous measurements and this extended emission 
is entirely contained within the 60 arcsec aperture used to calculate 
the photometry. Fluxes are in good agreement with the measure- 



ment of the peak flux b y ,SandelIl([2003[); [Siringo et al.[p009[ ) and 



the integrated flux by [Jenness et al. p0()2 1 



6.3 Arp220 

Arp220 is an ultra-luminous nearby galaxy and is well-studied 
in the submillimetre. Several publications have investigated the 
continuum flux at 850 |j,m {[Dunne et al. 2000 , L isenfeld et al.[ 
[20001 ISeiffert et ari[2507[ [Sakamoto et aL..2008| ). At 850 |xm, 
we calculate an integrated flux of 0.81 zb 0.07 Jy, compared with 
previous determinations ranging from 0.78 Jy (Sakam oto et al.[ 
(2008 )) to 0.83 Jy ( Dunne et al.K 2000)). At 450 |^m, the integrated 
flux of 5.4 Jy agrees reasonably with the SED fit from [Truch et aL[ 
( [2008| ). There are no published continuum fluxes for this source at 
450 t^m. [Truch et"aLl ( [2008] ) fit an SED to Arp220 using a combi- 
nation of published submillimetre and far-infrared flux densities, to 
calculate interpolated flux densities in their 250, 350, and 500 ^im 
BLAST bands (the last of which overlaps significantly with the 
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Table 4. Measured fluxes of SCUBA-2 secondary calibrators 



SOURCE 


RA 


Dec. 


N(obs) 


850[I] 


450[I] 


850[P] 


450[P] 


Note 




(J2000) 


(J2000) 




(Jy) 


(Jy) 


(Jy beam~^) 


(Jy beam~^) 




CRL618 


04:42:53.67 


+36:06:53.17 


113 


5.0 ± 0.2 


12.10 ± 1.05 


4.89 ± 0.24 


11.50 ± 1.4 


y 


CRL2688 


21:02:18.27 


+36:41:37.00 


173 


6.13 ±0.211 


29.10 ±2.5 


5.64 ± 0.27 


24.9 ± 2.9 


y,(a) 


Arp220 


15:34:57.27 


+23:30:10.48 


20 


0.81 ± 0.07 


5.4 ± 0.7 


0.79 ± 0.9 


5.2 ±0.8 


y 


PV Cep 


20:45:54.39 


+67:57:38.8 


21 


1.35 ±0.05 


10.6 ± 0.86 


0.82 ± 0.04 


5.7 ± 0.65 


n,(a) 


MWC 349 


20:32:45.53 


+40:39:36.6 


14 


2.19 ±0.08 


3.2 ± 0.25 


2.21 ±0.11 


3.4 ± 0.26 


n 


HD169142 


18:24:29.78 


-29:46:49.37 


12 


0.58 ± 0.02 


3.41 ± 0.24 


0.52 ± 0.03 


2.21 ± 0.25 


y, (a,b) 


V883 Ori 


05:38:19 


-07:02:2 


3 


2.00 ± 0.07 


10.4 ± 1.0 


1.55 ±0.09 


7.8 ± 1.00 


y, (b) 


HLTau 


04:31:38.44 


+ 18:13:57.65 


9 


2.42 ± 0.08 


10.3 ± 0.86 


2.32 ±0.1 


8.3 ± 1.03 


y, (b) 


BVP 1 


17:43:10.37 


-29:51:44.00 


4 


1.55 ±0.05 


16.9 ±0.8 


1.37 ±0.07 


11.5 ± 1.5 


y,(a,b) 



[I] indicates integrated flux measurements. [P] indicates peak flux measurements. Error is quadrature addition of uncertainty in the FCF 
measurement and the rms noise in the co-added map for each calibrator. Notes: y/n indicates if the source will continue to be observed as 
a SCUBA-2 calibrator, (a) Extended emission at 450 |Lim. Col. 4 gives the number of observations taken per source between 201 1 July and 2012 
May. (b) Low S/N measurement, particularly at 450 |j.m. Accuracy of these flux densities is uncertain, and further observations are required. 
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Figure 12. 850 |Ltm (left) and 450 |xm (right) image of Arp 220. The colour 
scale is linear as shown in the colour bar on the right hand side. In each map, 
the five contours scale linearly from the peak signal as shown in Table [4] 
down to the 3-cr detection level (0.18/0.36 Jy at 850/450 |xm respectively). 



Figure 13. 850 |j.m (left) and 450 |Lim (right) image of PV Cep. The colour 
scale is linear as shown in the colour bar on the right hand side. In each 
map, the five contours scale logarithmically from the peak signal as shown 
in Table|4]down to the 3-cr detection level (0.03/0.22 Jy at 850/450 |xm re- 
spectively). 



SCUBA-2 450 |^m bandpass). The resulting fluxes derived from 
20 observations of the source in 2011 June and July show good 
agreement with previous results. From a selection of observations 
in the best conditions, this unresolved source shows a measured 
FWHM of 13.1 arcsec at 850 |^m and 7.95 arc sec at 450 ^im. This 
agrees with the main-beam Gaussian fit in Table [T] Though faint, 
Arp 220 will continue to be observed as a calibrator for SCUBA-2 
in good weather. 



6.4 PV Cep 

PV Cep is an irregular variable star, with strong extended emission 
at both SCUBA-2 wavelengths. This is reflected in the increased 
integrated flux in comparison with its peak value. It is, for 
these reasons, unlikely to be used as a standard calibrator. The 
peak fluxes at both wavelengths are significantly lower than the 
background-corrected flux densities reported by [Sandell et al.| 
(j2011| and this is quite likely a result of the variability of the 
source. [Sandell et aL| ([2011) also postulate that their background 
fitting at 450 \im might produce elevated 450 \im results. Figure pj] 
shows the well-resolved image of the circumstellar disk, and 
surrounding extended emission. PVCep will not be a regularly 
observed calibrator as CRL 2688 offers a better alternative in this 
part of the sky. 
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Figure 14. 850 M-m (left) and 450 |Lim (right) image of MWC 349. The 
colour scale is linear as shown in the colour bar on the right hand side. 
In each map, the five contours scale linearly from the peak signal as shown 
in Table|4]down to the 3-cr detection level (0.03/0.18 Jy at 850/450 ^m re- 
spectively). 



6.5 MWC 349 

MWC 349, while quite close to CRL 2688, has the advantage 
of being a point source object. Rejected as being too faint for 
calibration with SCUBA, it is satisfactorily detected at both wave- 
lengths with SCUBA-2 as shown in Figure [14] and unresolved, 
with measured FWHM of 13.2 and 7.9 arc sec at 850/450 M,m 
respectively, agreeing well with the main-beam Gaussian fit in 
Table [1] At 850 \im, we measure a flux of 2.19 ± 0.08 Jy, lower 



© 2012 RAS, MNRAS 000,[l]{T2] 



10 J.T. Dempsey et al. 




Right Ascension Right Ascension Right Ascension Right Ascension 



Figure 15. 850 |xm (left) and 450 ^tm (right) image of HD 169142. The 
colour scale is linear as shown in the colour bar on the right hand side. In 
each map, the five contours scale linearly from the peak signal as shown in 
Table|4]down to the 3-o- detection level (0.03/0.17 Jy at 850/450 \im respec- 
tively). 



by a statistically significant amount to the 2.6 ± 0.07 Jy seen 
by [Sandell et al.| ( |2011| ). It is an emission-line star with strong 
free-free emission. We also detect a lower 450 \xm flux of 3.2 ± 
0.25 Jy, in comparison to the 5.0 zb 1.1 Jy SCUBA observation 
reported by [Sandell et aL|p011| ). [Sandell et al.| ( |2011| ) suggest that 
the background- subtraction method adopted with the SCUBA data 
produced questionably high 450 \xm fluxes in some cases, though 
this does not explain the flux discrepancy at 850 |^m. As these 
data were taken in a short period and with a long interval since 
the SCUBA measurements, source variability cannot be ruled out. 
CRL2688 remains a preferred calibrator to this source, which will 
not be regularly observed. 



6.6 HD 169142 

HD 169142 is a well-studied Herbig AeBe Star, and was selected 
as a study source for SCUBA-2 calibration to probe the sensitivity 
limits of the instrument, particularly at 450 \im. [Sandell et al.[ 
pOTT] ! find that the disk may be resolved, and report a flux density 
of 0.57 Jy at 850 ^tm and 3.34 Jy at 450 |im. We measure an 
asymmetric disk with a fitted FWHM of 15.5x13.4 arcsec at 
850 \xm which agrees with this finding. This is the faintest source 
observed with SCUBA-2 in this study and the 850 \im flux of 
0.58 Jy and 3.41 Jy agrees with the value measured by [Sandell] 
jet al.[ ( [201 1| ) within the uncertainty. The co- added maps of the 
source are shown in Figure p3] 



6.7 V883 0ri 

V883 Ori is FU Orionis star, which has a similar right ascension to 
a preferred calibrator, CRL618 but has the advantage of being in 
the south. [Sandell & Weintraub[ ( |200Tl ) observed this with SCUBA 
and found it to be unresolved, with little evident extended emission 
(though they suggest a possible emission region to the eastern 
edge of the 850 |^m map). They measured a background-corrected 
flux density of 1.41 ± 0.02 Jy at 850 |^m and 9.45 ± 0.17 Jy 
at 450 \im. Additional observations with SCUBA measured an 
average peak flux at 450 \im flux of 7.28 zb 0.07 Jy ( |Barnard[ 
[2005] ). We measure an 850 |^m integrated flux of 2.0 zb 0.07 Jy 
and a 450 ^im integrated flux of 10.4 ± 1 Jy. [Siringo et al.[ ( [2009| ) 
reports an integrated flux from LABOCA at 870 |j,m of 1.93 Jy, 
which agrees within our error range. The SCUBA-2 and LABOCA 



Figure 16. 850 |Ltm (left) and 450 |Ltm (right) image of V883 Ori. The colour 
scale is linear as shown in the colour bar on the right hand side. In each map, 
the five contours scale linearly from the peak signal as shown in Table |4] 
down to the 3-c7 detection level (0.14/0.37 Jy at 850/450 jiim respectively). 




42.0 40.8 39.4 4:31:38 37.20 36.0 34.80 42.0 40.8 39.4 4:31:38 37.20 36.0 34.8 
Right Ascension Right Ascension 



Figure 17. 850 |j.m (left) and 450 |Lim (right) image of HL Tau. The colour 
scale is linear as shown in the colour bar on the right hand side. In each map, 
the five contours scale linearly from the peak signal as shown in Table |4] 
down to the 3-cr detection level (0.08/0.44 Jy at 850/450 |Ltm respectively). 



measurements both disagree with the lower 850 ^im integrated flux 
measured by [Sanden| ( |2003| ) by a statistically significant margin. 



6.8 HLTau 

HL Tau is a young T Tauri star, and was used as a secondary 
calibrator for SCUBA fS andell[|2003l jJenness et aLl|2002l ). It is 
well studied in the submillimetre, particularly as imaging suggests 
a low-mass companion object in the parent disc material ( [Greaves j 
jet aL|2008| ). The SCUBA-2 maps here do not detect a significant 
amount of this disc with a measured FWHM at 850 \im of 13.8 arc 
sec and 8.3 arc sec at 450 \im. The integrated fluxes agree well at 
both wavelengths with the measurements of jJenness et al.[ ( [2002| l. 
Only a small sample of data was taken of this source in the period 
available. Further observations will be conducted to explore the 
properties of the extended emission with better S/N. 



6.9 BVP 1 

Figure[T8]shows the SCUBA-2 images of BVP 1, a bright, compact 
submillimetre source located near the galactic centre. SCUBA 
observations of the source are to date the only published results 
at these wavelengths, and II Barnard et al.[ p004| ) measured a peak 
flux at 1.069/1 1.318 Jy at 850/450 ^im respectively. It is a good 
potential secondary calibrator for SCUBA-2 as there is a paucity 
of good calibration candidates in this region of the sky. With only 
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Figure 18. 850 |xm (left) and 450 |xm (right) image of galactic source 
BVPl. The colour scale is linear as shown in the colour bar on the 
right hand side. In each map, the five contours scale logarithmically from 
the peak signal as shown in Table |4] down to the 3-cr detection level 
(0.144/0.46 Jy at 850/450 |xm respectively). The correlated background 
structures are likely a result of the proximity of the source to the galactic 
centre. 

a handful of observations to date, SCUBA-2 measures peak fluxes 
of 1.37/11.9 Jy at 850/450 ^im respectively. This is higher at both 
wavelengths than Barnard et al. (2004), with a 4-cr discrepancy 
at 850 M^m. The integrated fluxes are 1.55/16.9 Jy at 850/450 M,m. 
Fits of the source profile show that it is compact at 850 |^m, with 
a central FWHM of 14.3 arcsec. It shows extended emission at 
450 i^m, as indicated by the elevated ratio of integrated to peak 
flux, though its central core is extremely bright and fairly compact 
with a measured FWHM of 8.8 arcsec. The SCUBA observations 
were taken in jiggle-map mode, and in the crowded region near 
the galactic centre, with large amounts of extended emission, it is 
possible that chopped emission resulted in the lower flux recorded 
by [Barnard et al.| ( |2004| ). Observations and investigation of this 
object continue and will hopefully shed more light its properties. 



7 CONCLUSIONS 

SCUBA-2 is fully commissioned and currently running very suc- 
cessful science operations on the JCMT. Atmospheric calibration 
at both wavelengths has been optimised using new line-of- sight 
opacity algorithms and the derived atmospheric extinction relations 
are shown here. Using a large sample of primary and secondary 
calibrator observations, the flux conversion factors (FCF) for the 
850 i^m and 450 ^im arrays are deduced and an investigation of 
the instrument beam- shape and photometry methods is given. 
The relative flux calibration is better than 10 per cent at both 
wavelengths. The sample-size of the calibration observations and 
accurate FCFs have allowed the determination of the 850 ^tm and 
450 \xm fluxes of several well-known submillimetre sources and 
these results are presented. 
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